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Abstract 

We present a much simplified version of the embedded point mass lens theory and then make the 
obvious extension of this theory to any embedded transparent lens. Embedding a lens effectively 
reduces the gravitational potential's range, i.e., partially shields the lensing potential because 
the lens mass is made a contributor to the mean mass density of the universe and not simply 
superimposed upon it. This presentation results from simplifying previous formulations of the 
embedded point mass theory and is quite similar to the standard theory. We are then able to re- 
derive the simplified point mass theory by applying Fermat's least-time principle to the time-delay 
function. Even though rigorous derivations are only made for the point mass, the lens equation 
to the lowest order for any distributed lens, is obvious. We find from this simplified theory that 
embedding can introduce corrections at the few percent level in weak lensing shears caused by 
large clusters but only at large impacts. The potential part of the time delay is also affected in 
strong lensing at the few percent level. Additionally we are surprised to find that the cosmological 
constant is only present in this theory at a higher order. 
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In this Letter we present a simplified version of the embedded point mass lens theory jl-4| 
that can easily be used by anyone familiar with conventional lensing theory, see Eqs. (15])- 
(fTUj) . Because we have succeeded in replacing image impact position variables by the image's 
angular position Qj (Fig. 1), we have succeeded in giving the point mass lens equation in a 
form which is easily compared with the conventional theory [see Eq. ([13]) ]. Image properties 
such as the time delay, i.e., the excess time light takes to reach the observer because it 
encounters the embedded lens, are also given. By applying Fermat's least time principle 
5|, |6[ to this time-delay function the lens equation to lowest order is then consistently arrived 
at [see Eq. i.e., it is the same as the lowest order term in the rigorously derived result, 
Eq. (I13p . Because we can interpret this lowest order equation within the conventional theory 
as lensing by a point mass plus a negative surface mass density [see Eqs. (Q and (I14p ]. we 
are able to extend this theory to any embedded transparent lens, see Eqs. (I15p and (fLSj) . 
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In previous work jl-4| we have investigated, in depth, effects of embedding on the point 
mass lens equation and on all resulting image properties. However, application of that 
theory is quite complex because of the lens impact variables used. Our embedded lens 
model remains a spherical comoving mass void (see Fig. 1) in a flat pressureless Friedman- 
Lemaitre- Robertson- Walker (FLRW) universe, containing the cosmic mass m removed from 
the void condensed at its center [?], lsj]. The actual geometry in the void is described by the 



Kottler metric |s| which is similar to Schwarzschild's geometry but with the addition of a 
cosmological constant A. This lens model is used to maintain absolute correctness of the 
gravity theory and hence avoid approximation errors in the lensing theory. A universe filled 
with such voids is called a Swiss cheese universe and has historically been used to compute 



10Ml5| and on the anisotropy of the 



effects of inhomogeneities on distance-redshift relations 
CMB 16). The comoving radius Xb of a lens void at redshift Zd is related to the Schwarzschild 
radius r s = 2Gm/c 2 of the embedded lens by 



n^xl (1) 



where Hq is the Hubble constant, Q m is the usual matter density parameter, and the current 
radius of the flat universe has been taken equal to 1. The lens void expands with the 
background cosmology and had a physical radius r d = Xb/i^ + Zd) when the observed photons 
passed the deflector. Throughout a subscript 'd' on evolving quantities means evaluated at 
the deflector's redshift z^. We start by presenting the new point mass results, explain how 
they were arrived at from the previous more complicated expressions, and then explain how 
they are generalized to distributed lenses. 

Only two parameters appear in the simplified point mass theory, the familiar angular 
Einstein ring radius 

and the maximum value of the angular image position dj (see Fig. 1) 

0m = ^ = ^- (3) 
D d Xd 

where Xd is the FLRW comoving distance to the deflector from the observer and Dd, D s , 
and D ds are angular diameter distances, respectively of the deflector, source, and source 
relative to the deflector. The maximum image angle 9m occurs for the primary image when 
the light ray from the source just grazes the spherical void and can be thought of as the 
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FIG. 1. A photon travels from a source located at angular position 0$, a comoving distance Xs 
from the observer, and enters a Kottler void of comoving radius Xb centered at comoving distance 
Xd from the observer. The photon is deflected by an angle a (< 0) and returns to the FLRW dust 
on its way to the observer, where it appears at angle 9i. The maximum image angle Q\ is 9 m and 
occurs when the light ray grazes the Kottler void. 



angular size of the embedded void as seen by the observer. For this maximum image angle, 
the deflection angle is zero and the source is likewise located at 9m- For a massive galaxy or 
a cluster lens (M = 10 12 M Q or 10 15 M Q ) at redshift Zd = 0.5, and a source at z s = 1.0, the 
Einstein ring radius 9e is respectively 1".66 or 52'.'6, much smaller than the maximum image 
angle, i.e., 0e/0m — 0.0085 or 0.027. Consequently both strong and weak lensing effects 
are seen within the void. By using 9j as the independent image variable we find from our 
previous work fl-fl that up to a constant, the time-delay function cT (accurate to order 8 , 
explained below) is 



cT(9 s ,9j 
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This expression, without the (1 + Zd) factor, is often called the Fermat potential, see (17J. It 
represents the difference in arrival times of two signals originating at a source at redshift z s , 
one traveling on a straight path entirely in FLRW and the other encountering the embedded 
deflector at redshift z d (see Fig. 1). It is the sum of a geometric and a potential term 
written in a form to which Fermat's principle can be applied 5|, |6fl. By varying cT(9s,9i) 
with respect to the image angle 9j, i.e., by applying Fermat's least time principle, we obtain 



the embedded lens equation, accurate to order S l , 

9 2 



E 



9 s -0i- 9j 



v / HW • (5) 



The conventional linear lensing result, without embedding, is regained by neglecting the 
(0i/6 M ) term. 

From Eq. (jSJ) the inverse image matrix eigenvalues (respectively axial and radial) are 
found to be 

0s -, f 0e 



0i ' 

(0 ' 



(^) [V^-(0i/0m) 2 ]\ (6) 

r t = i + (^) 2 [1 + 2 (0i /0m) 2 ] a/i - (0i/0m) 2 - (7) 



dd 

The reciprocal amplification is 
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lT x = H a r = 1 + 3 (J^) Vl-^/M 2 " (^y) [1 + 2(^ A/ ) 2 ] [1 - (^/^/) 2 ] 2 



with a normalized effective surface mass density k present 
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K = 1 - ( H + a r )/2 = -- i-^-j y/1 - {erle M y. (9) 

This negative term can be interpreted as the void mass removed from the observing light 
beam at image position 6j. The shear in the bundle of light rays at image position 9i is 

7 = K - a*)/2 = \ {j-^j [2 + {0i/0m?WI-{0i/0mY. (10) 

By neglecting the 0e/0m and QjjQu terms, i.e., the shielding terms caused by embedding, 
the above expressions all reduce to the conventional non-embedded values. These results 
represent a significant simplification in our previous work on embedded lenses and can 
easily be used without following the complexities of their derivation outlined next. The 
generalization of the above point mass lens equation, Eq. (JSJ), to spherically distributed and 
arbitrarily distributed lenses is given respectively in Eqs. (I15p and ( TT6|) . In Fig. 2 on the left 
we have plotted the shear 7 as a function of 8j and have compared the embedded 7 with the 
conventional value for a large cluster lens. The shear is seen to differ from the conventional 
Schwarzschild value by more than a few percent but only at impact angles more than 15 @e- 
In Fig. 2 on the right we have plotted time delay differences caused by embedding. These 
differences are primarily due to the potential part of Eq. (j4j) differing from the Schwarzschild 
value at large primary image angles. 
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FIG. 2. (Left) Embedding corrections to the image shear 7 for a 10 15 Mq cluster lens. The source 
and lens redshifts are z s = 1.0 and = 0.5, respectively. The solid blue curve is the shear, 7, and 
the dashed red curve is the fractional difference in the image shear, A7/7, caused by embedding. 
Both are plotted as a functions of 6j/6m where the Einstein ring radius is Be =52'.'6, the void 
radius is 9 m =32'.8, and 9e/9m = 0.027 is the vertical dotted line. (Right) The difference in 
the embedded time delay and the conventional Schwarzschild time delay, AT — ATs c h, computed 
using Eq. (j3|) and the Schwarzschild equivalent for the primary and secondary images, is plotted 
as a function of source position angle 9s /6m- The time delay difference is given in months for 
four point mass lenses, m = 1O 15 M , 1O 14 M , 1O 13 M , and 10 12 M Q (respectively top to bottom 
curves). The small vertical ticks on each curve are at respective 6e/9m values. 



The rigorously derived embedded point mass lens equation, Eq. (1131) . and from it the 



2J where we computed the lens 
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lowest order term Eq. ([5]), follows from our earlier work 

equation and other lensing quantities as functions of two impact variables of the 
metric, the minimum impact radius tq and the void impact angle 0i (see Fig. 1 of [4|). In 
3|, |4| we succeeded in eliminating the minimum impact radius r in favor of the single void 
impact angle 0i and were consequently able to write the lens equation, etc., as functions of 
4>i alone. Only now have we succeeded in eliminating <pi in favor of observable 61 and truly 
simplifying the point mass embedded lens equation as given in Eq. (fT3|) . 

The conventional linear lensing theory is only an approximate theory for a variety of 
reasons including the fact that its two fundamental ingredients, the bending angle a = 
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-2r«/ro, and the image position 9j = ro/D d , are both lowest order approximations, see 
2j. Obvious corrections to a require terms of order (r s /ro) 2 and geometry corrections are 
required for 9j. In addition to these obvious corrections we have investigated corrections that 
come from insisting that the gravity theory confirm to Einstein's. Such insistence requires the 
deflector contribute to the mean cosmological mass density, i.e., that it is embedded in the 
background cosmology. To track all corrections we have ordered them by their magnitudes 
using a parameter 5, whose value is 1. In jl-4] the details of tracing a ray from a source at 
redshift z s through its interaction with an expanding void at redshift z d was done using a 
series of approximation of orders up to d 4 to obtain the needed accuracy, e.g., to evaluate 
;he time delay as given in Eq. (EE)), see Eqs. (10) and (14) of 2| and Eqs. (22) and (23) of 
4|. When observing galaxies and clusters at cosmological distances, lensing parameters of 
order S 1 are the angular void radius 9m, see Eq. (j3J), and f3 d , the expansion rate of the void's 
boundary as measured by static observers in the void 



d = ^7 = \jr s /r d + Ar 2 /3 = ^r d , 



(11) 



where H d is the Hubble parameter at the lens redshift z d . Equation (flip , as well as Eq. (TfJ), 
are consequences of the geometry of embedding. When r s /r d and appear individually 

they were both treated as terms of order S 2 . For the galaxy and cluster lens discussed above, 
S 1 terms are of respective orders ^ 10~ 3 and < 10 -2 . 



The steps we fol 
inverting Eq. (8) of 
9i, 



owed to obtain the simplified embedded theory were as follows: By 
4J we obtained the previously used impact variable <pi as a function of 



H = sin- 1 {9 I /9 M )+S/3 d (9 I /9 M )S 



2 0\i (Oi/OmY 



+ 



a/1 - (0i/0m) 



-0(8 3 ). 
(12) 

The embedded lens equation was then obtained by eliminating (f>i from Eqs. (6) and (7) of 
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4J using Eq. flTJ 
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The simplified lens equation ([S]) is obtained by dropping all correction terms of order 8 2 
and higher in Eq. ([TBI . To the order of accuracy needed to obtain Eq. (j3J) it is sufficient 
to replace the angle 0i used in {4] by 9j using only terms up to order S 1 of Eq. (1121) . The 
geometrical part of the time delay given in Eq. fll]) was obtained by using Eq. (10) of ^ 
combined with Eq. (7) of 4j. The potential part of the time delay given in Eq. (Ill) was 
obtained by eliminating 0i in Eq. (23) of jjj. The simplified embedded lens equation ([5]) is 
then re-affirmed, as expected, by applying Fermat's least time principle to Eq. 01]). Knowing 
Eq. (|5]) is sufficient to reproduce the image properties for the embedded point mass accurate 
to order S 1 . 

The absence of terms of order S 1 in Eqs. (fT3|) and (jlj) was surprising since they were 
previously present when other impact variables were used. When the lens equation and the 
time delay were previously given in l|-|4| as functions of tq and/or <j>i, (3d terms (terms of 
order S 1 ) were actually present just as they were in the bending angle a. What we see in 
Eqs. (Hj)- (|10p is that when 9i is used as the independent variable, terms of order S 1 all vanish, 
including fa terms, and hence A is no longer present at this order. We find this somewhat 



surprising in light of an earlier controversy 



19H22J, settled in the affirmative, about the 



existence of an effect on the bending angle a by the cosmological constant. However, we do 
find that when terms of order S 2 are kept in the lens equation (parameterized by 9j), A is 
once again explicitly present |23|, [24|, see Eq. (fl3~|) . 

The point mass theory is of limited usefulness on the scale of clusters. Fortunately the 
above is easily generalized to lenses with arbitrarily distributed mass, once it is observed 
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that the order S 1 lens equation (jSJ) is the same as the classical lens equation for the sum of 
a point mass and a distributed mass lens with negative projected surface mass density 

c 2 D, 



Yj rr K 



-K. 



(14) 



AuGD d D ds 

Here k is from Eq. ([9]) and S cr is the conventional critical surface mass density [17|. This 
leads to the immediate observation that we can generalize, to this order of accuracy, lens 
equation (J5J) to include any transparent spherically symmetric mass 

n 3 N 
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where f(9j) is the fraction of the actual distributed mass lens within the impact cylinder 
defined by_ Oj. Because an exact Einstein model is known, the Lemaitre-Tolman-Bondi 



metrics 25 



- 27] embedded in FLRW, this result can be confirmed by rigorously extending 
the work done in [l-4]. The generalization to any distributed lens, not necessarily spherically 



symmetric, is also obvious 

D ds 
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where ct is the conventional bending angle without embedding. This result cannot be readily 
confirmed because the relevant Einstein solution is not known. Nonetheless, it should be 
correct to order S 1 . Image properties for these transparent lenses can be obtained much as 
they were for the point mass in Eqs. fl6|)- (|T0|) . Because of the simple dependence of these 
lowest order embedded lens equations on angular size distances, they should also be valid 
for any FLRW background models, not just flat cosmologies. To extend the higher order 
lens equation (TT5j) to non-flat backgrounds will require additional work. 
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